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1 Second order ordinary differential equations 



invariant under time-translation and rescal- 



ing 

The general form of the second order ordinary differential equation invariant 
under the two symmetries representing invariance under time translation and 
rescaling, viz 

is 

where / is an arbitrary function of its argument. As a second order equation 
with two symmetries (|1.2| ) can always be reduced to quadratures. In the 
particular case that / is linear with a specific form, ie 

/ = 36(|)+6^ (1.3) 

giving the equation 

x + Sbxx + b'^x^ = (1.4) 

the equation possesses eight Lie point symmetries, naturally with the algebra 
sl{3,R), and so is linearisable |1^. 



Not only does ( |1.4|) possess the Painleve property but it has both Left and 
Right Painleve Series ||1T|. The equation is a representative of the Riccati 
hierarchy 0|[6.39, p 550] and is integrable via transformation to a linear third 
order equation using a generalised Riccati transformation. We write 

X = a—, (1.5) 

w 

where a is a constant to be determined. Equation (|1.4| ) becomes 

W , , , WW /,9 9 , \ w^ , ^ 

— + 3 a6 - 1 W + ih^a^ - 36a + 2 ) — = 0. 1.6 

w w^ ^ ^ 
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Equation ( |1.6| ) is a third order hnear differential equation in the case that 
both 

a6 - 1 = and fe^a^ - 36a + 2 = 0, (1.7) 
ie ab = 1. The solution is 

w = A + 2Bt + Ct'^, (1.8) 
where A, B and C are arbitrary constants. Consequently 

x(t) = -- 7 , , 1-9 

where ki and /c2 are arbitrary constants. The second of (|1.7| ) is zero if ah = 2. 
In this case (|L^) becomes 

ww + 3ww = (1-10) 

which is a linear third order equation in the variable w"^ with the same solution 
as ( |T78|) with the exception that the w is replaced by w^. In either case the 
solution of ( |1.4| ) is the same when the transformation ( |1 . 5| ) is applied. It is 
amusing to note that, if (|1.10|) is multiplied by the integrating factor w'^ and 



integrated, we obtain 

w = ^, (1.11) 

where K is the constant of integration, which is a particular instance of the 
Ermakov-Pinney equation Thus one can easily see the relationship 



between the solutions of the two third order equations. The solution of the 
Ermakov-Pinney equation 

X + Uj'^X=ir (1.12) 
X"^ 

is given by 

X = ^Au^ + 2Buv + Cf2, AC-B'^ = W, (1.13) 
where u and v are two linearly independent solutions of 

x + cj2a; = (1.14) 

and W is the constant value of their Wronskian. Without the constraint on 
the constants of integration the expression in the square root is simply the 
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general solution of a third order equation of maximal symmetry, the solution 
of which can be expressed in terms of the three independent quadratic terms 
obtainable from the solution set of the corresponding second order equation 



We can transform the second order differential equation ( |1.4|) into a two- 
dimensional system of first order equations by a transformation which pre- 
serves the symmetries in We write 




:i.i5) 



to obtain the system 



X = ax + /3xy (1-16) 
y = ~{a + b){2a + b)x^ -3{a + b)xy - f3y'^ (1.17) 

in which the constants a and f3 are arbitrary. 

The generality of ( |1.4| ), ( |1.16|) and ( |1.17|) is merely apparent. In ( |1.4] ) 
the parameter b may be scaled out of the equation by the transformation 
bx ^ X or, alternatively, by a rescaling of the time variable. In the case of 
equations (|1.16|) and ( p..l7|) we can remove the parameter f3 by the rescaling 
of the new variable y according to the transformation (3y — > y. Consequently 
the analysis of these equations can be performed without loss of generality 
on the second order ordinary differential equation 

x + 3xx + x^ = (1.18) 

and the system of two first order equations 

X = ax^ + xy (1-19) 
ij = -(a + l){2a+l)x^ -3{a + l)xy -y^ (1.20) 

obtained from ( |1.18D by the transformation 

y = ^-ax (1.21) 
x^ 

in which we have replaced the parameter a by the parameter a to aviod 
confusion when we use the former in the Painleve analysis. Note that G2 of 
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the original second order equation becomes 

Odd 

for the two-dimensional system. 

The two-dimensional system, (|1.19| ) and ( |1.20|) , is a particular instance of 



a class of quadratic systems integrable by the factorisation method of Adler- 
Kostant-Symes and its generalisation P, ^] which have been recently studied 
from the point of view of their symmetry and singularity properties [jlO . 



2 Explicit integration of second order ordi- 
nary differential equations of the class 

We consider certain specific forms of (|1.2|) , viz 



4 + /f4)=0, x = xit). (2.1) 



2.1 A) f{x/x^) = 2ax/x'^, a real. 
Equation (|0| ) becomes 



X + 2axx = (2.2) 
and, since ( p. 2]) can be written as 

dp 

— = —2ax, pix) = X, (2.3) 
ax 

by carrying out the trivial integration of (p.3[) we finally obtain for x{t) 



Equation (^I^) has a simple pole as singularity. 
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2.2 B) f{x/x^) = 2ax/x^ + 26, a and b real. 
In this case ( pTip becomes 

X + 2axx + 26x^ = 0. 



(2.5) 



Equation ( [^.5D appears in Kamke |0|[6.43, p 551], but the method used 
there leads, even for special values of a and b, to integrals which are beyond 
any hope of calculation. Therefore we resort to transforming (|2.5|) into a 



third-order equation by means of the Riccati Ansatz |11 



w 



X = a—, w = w{t), 
w 



in which case ( p.5|) becomes 



(2.6) 



W WW 



w 



+ ^(-3 + 2aa) + —(2 - 2aa + 2ba 



w w^ 



w" 



(2.7) 



Choosing in (|2.71) 



a 



2a^ 
9 



(2.8) 



we finally obtain from ( p.7| ) and ( |2.6| ) that the general solution of ( |2.5| ) is 



x(t) = - 



provided 



2a' 
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(2.9) 



(2.10) 



where ki and /c2 are constants of integration. Equation (|2.9| ) has poles as 
singularities. The results ( [^.9| ) and ( [^.10 ) for ( p^.5| ) are known |TT|| and are 
identical to ( |1.9| ) for (|1.4|) provided that we identify the a in (|2.5| ) with 36/2 
in (0). 

In the following we drop the condition ( |2.1CI| ) and investigate the resulting 
third-order equation (p.7|), ie 



9b 

w'^ w + w'^c = 0, c = — - — 1. 

2a^ 



(2.11) 
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We reduce the order of by using the standard substitution 



w = p, p = p{w). 
Hence, from ( ^.11| ) and ( |2.12| ), we deduce that 

2 // I 2 /2 I 2 n 

w pp + w p + cp = (J 



(2.12) 



(2.13) 



which is a hnear second order ordinary differential equation in the variable 
p^ of Euler type. Consequently ( p.l3| ) has eight Lie point symmetries which 
is an increase from the two of the original second order ordinary differential 
equation. The technique of increasing and decreasing the order of an equation 
using different symmetries has been found often to be most useful |jl|, |^ . 
Equation (|2.13|) is treated in Kamke P] [6.184, p 585] and the result is 



Pi = l{i- VT^) , p2 = |(i + v^r^' 

Equations (^J^) and {^J§ yield 



(2.14) 



c^l/8 



w 



-Pi/2 



Ci + C2W 



P2-P1 



-1/2 



dw = t + C, 



(2.15) 



where Ci, C2 and C are constants of integration. By means of the substitu- 
tion 

P2 - Pi n 5- 37 



w 



P2-P1 



-, yr^8c 



l-pi/2' " A--i 

the integral, /, on the left hand side of (|2.15| ) becomes 

2 



2-pi 



(2.16) 



(2.17) 



Now, / in ( p.l7| ) can be transformed into an integral of rational functions 
[|][pp 70 - 71] if 1/7 = integer = k. The case k = 1 leads, due to ( p.l4|) and 
( |2.16|) , to \/l — 8c = 1, c = 0, which in conjunction with ( p.lO|) and ( p.llj ) 



(2.18) 



produces the already known result, (|2.9| ). Let then k ^ 1 and set 
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By means of ( |2.18| ) the integral in ( p. 171 ) is written as 



C27(2 - Pi) 



fe-i 



du. 



(2.19) 



In ( p.l9| ) k — 1 = integer. By virtue of ( p.l6| ) < 7 < 4. Thus from 1/7 = 
integer = k we have k = 1,2,3, .. .. Let k = 2, ie = 1/2 . Then by ( p.llj ) 
and ( CT) 



and from ( p.l9| ), on inserting pi from ( p.l4| ) and ( |2.16| ) with 7 = 1/2, and 
utihsing ( p.l5|) , we obtain 



- 3Ciu 



9(t + C)C| 
14 



0, u = u{t). 



Owing to (pl6D, ( |2.16| ) and ( p.l8[ ) we obtain, since 7 = 1/2, 

7uu 



x{t) 



a (m2 - Ci) ' 



(2.21) 



(2.22) 



where u{t) is given imphcitly by ( p.21| ). 

Equation (|2.21| ) can be solved analytically for u = u{t). To this end we 
introduce the notations 101 



-cl + ^,ct{t + cf 



1/2 ^ 1/3 



(2.23) 



and 



S2 
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28 



Cl{t + C) 



-cl + ^,C',{t + Cf 



1/2- 



1/3 



(2.24) 



Then, from (|2]2T|) 



M = 5*1 + 5*2 

^ = -i (51 + ^2) + ^ (51 - 52) 

« = -i(5i + 52)-^(5i-52) 



(2.25) 
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and, after a modicum of calculation, from ( p.22| ) - ( |2.25| ) we deduce, consid- 
ering firstly M = 5*1 + 5*2 in (|2.25|) , 

x{t) = ^^^^^ (2.26) 

3a{ul - 1)^-1 + P(t + C)2 



where /c=^C'2/^i^^ ^^'^ ^ constants of integration, and 



s 

S2 



, = {k{t + C) + [-1 + k\t + Cf]^'^}'"" 

{k{t + C)-\-l + k\t + Cn''}"' (2-27) 

Ui = Si + S2 



Y{t) = k{sl~sl). (2.28) 

It can be easily verified that uf ^ 1. Consequently, the exact solution, 
( |2.26|) and ( p. 271) , of (|2.5| ), proviso (|2.20|) , possesses branch points as singular- 



ities. We have restricted ourselves to one of the three possible solutions u{t) 
of (|2.21|) appearing in ( p.25|) in order to gain some insight into the structure of 



the solutions of ( p.5| ) without unduly complicating the respective calculations, 
provided (|2.2CI|) holds. This clearly implies that in the initial value problem 
consisting of ( p.5| ) endowed with the initial conditions xito) = xo,i(to) = ^^o, 
we take into account initial conditions such that the unique solution x{t), 
deriving of course from the basic equation ( p.21|) , can also be extracted from 
u = Si + 32- The remaining possibilities for u in ( |2.25|) will most probably 
lead, owing to the form of Si and 5*2, to x(t) possessing branch points or 
essential singularities or both. 

From (|2.19| ) it is obvious that for k = 3,4, . . . the x{t), which now holds 



for all a and b fulfilling ^T^Sc = 3/(4A; - 1), c = (96 - 2a^)/{2a^), is still 
being given by ( 2.22|) while u = u{t) is determined implicitly through an 



algebraic equation of degree 2k — 1, the solution of which for u = u(t) cannot 
be effected analytically. Nevertheless the aforementioned x{t) and u = u{t) 
do represent an exact solution of ( |2.5| ). The remaining case for which / in 
(|2.17|) can be transformed into an integral of rational functions |§][pp 70 



71] is 1/7 -1/2 = integer = k. If = 0, then 7 = 2 and by (|1|) 6 = 0, 
which is simply case A) above. As < 7 < 4, we have k = 1,2,3, .. . and we 
set 

Ci + C2V^ = uV, (2.29) 
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so that (g^D - (|2T7| ) yield 



7(2 -Pi) 



-(fc+i) 



du = t + C, 



(2.30) 



while the equivalent to ( p.22| ) is now 

x{t) 



3uu 



a(C2 



(2.31) 



Clearly the integral in (|2.30|) is in principle calculable, but even for k = 1 
the resulting defining equation for u = u{t) is not invertible. However, we do 
have in implicit form, ie ( p.30| ) and (|2.31|) , an exact solution to (|2.5|) provided 
VI - 8c = 3/(4/c + l), c = (96-2a2)/(2a2). We do not think it is worthwhile 
to pursue this further, since the characteristic features of the structure of the 
solutions to ( |2.5| ) are borne out by the cases analytically calculated above. 
Note also that, either by setting one of the two constants Ci and C2 in ( p.l5| ) 
equal to zero or directly from (|2.5| ), one immediately obtains the particular 
solution to (I2I 



VIZ 



x{t) 



A 



hy -a\+l = 0, 



(2.32) 



with one of the A- values in ( p.32|) for C = corresponding to (pl9| ) and ( p.lO| ) 
for ki = k2 = 0. The solution ( 2.32| ) has a simple pole as singularity. 

As can be expected, nothing is gained by interchanging pi and p2 in ( |2.14|) , 
while c = 1/8 in ( [^.14D leads to the appearance of logw in the integral in 
( |2.15| ) thus rendering it not calculable. This coincides with the result found 
when the Painleve test is applied to (|2.5| ) |1TT| . 

We conclude the discussion of ( |2.5| ) by noting that for a = 1/2 and 
b = —1/2 (|2^ ) possesses the general solution P][6.32, p 548] 



x{t) = Ci 



ViCit + C2;0,l) 
P(Cit + C2;0,l) 



V{Cit + C2; 0, 1) being the Weierstrass elliptic function [^[p 917], where the 
prime denotes differentiation with respect to Cit + C2, and Ci and C2 are 
constants of integration. By using the series expansion of V{Cit + C2; 0, 1) 
it is seen that x{t) in ( p.32| ^) has a simple pole at Cit + C2 = 0. 
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Finally, if in ( |2.5| ) we set 



X = u 



ISD becomes 



u - fw^u-i + 2auu~^/^ - 46 = 0. ( ^^ ) 

Then due to 

it = p, p = p{u) ( p.32[ ^) 

we deduce from ( p.32f ) 

upp - |p2 + 2apu^/'^ - Abu = 0. ([O^ 



Equation ( p.32p ) is an Abel's equation of the second kind and can be treated 
according to Kamke 0[pp 25-28]. The final result obtained after some cal- 
culations is 

(-2a)-i/2 / + Caz'-^-'-O"'^' dz = t + C3 

u = -2a {Ciz ^^^ + C2Z''')~\ z = z{t) (|23^ 
ri = (1 - ^1 - Ah/ a') /2, = (l + y'l -Ah/ a') /2. 

There are two independent constants in (|2.32f ) due to (|2.32p ). Equations 
( |2.32p ) are essentially equivalent to (p.l6|) - ( p.l9| ) and one obtains the same 



results as previously. Therefore transformation ( p.32f ) does not lead to any 



new results other than indicating another way to obtain the results already 
known. If instead of ( |2.32p ) we consider x = — equivalent results are 
extracted. 
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2.3 C) f{x/x^) = a{x/x^Y + hx/x^ + c, with a, h and c real 



and a 7^ 0. 

Now becomes 

XX + a(i;)^ + hx'^x + cx"^ = (2.33) 

Equation ( p.33| ) appears in Kamke [6.130, p 574] and we employ the 
method described there at the end of this paragraph. In the following by 
means of ( |2.6| ) we obtain from ( |2.33| ) by choosing 

a = — ^ (2.34) 



•■2 • S 

w w 
w + a h a— = 

w w\ 

d = a + 2 + c{ — - — j -(2a + 3). 



(2.35) 



We firstly assume 

c/ = (2.36) 

which implies one relation between the coefficients a, h and c. Then ( p.35| ) 
gives 

w = y 
y + a — = 

y 

the second of which is linear in 

Equation ( ^3^ ) is treated in Kamke [6.128, p 574]. We obtain, by 
taking into account (p.6|) after some straightforward calculation, 

2a + 3\ /a + 2 



fc2(t + fci)-V('^+i) + (t + fci) ' ^ ^ ^ 

as the solution to ( [^.33| ) if ( |2.36| ) is valid, ki and k2 being constants of inte- 
gration. If l/(a + 1) = -A;, A; = 1, 2, . . . or l/(a + 1) = fc. A; = 1, 2, . . ., the 
x{t) in ( p.38| ) possesses poles of order k at the most or poles of order + 1 at 
the most as singularities respectively. For all other values of l/(a + 1) ( p. 3^ ) 
may possess all kinds of singularities, for example a pole and an essential 
singularity if —1 — l/(a + 1) 7^ integer > 0, an essential singularity and a 
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branch point if l/(a + 1) + 1 7^ integer > and so on. In the case a = —1 
(|2.38|) becomes 

= F(It!*o' " = <'-^^' 

which has a simple pole as singularity. 

In what follows we discard condition ( p.36|) and examine the full equation 
(CT) - Then 

lb = p,p = p{w) 
w'^pp" + w'^p"^{a + 1) + dp^ = 

which is an equation of Euler type linear in the variable w°'~^'^ provided 
a 7^ —2. Equation ( p.40| ) is precisely of the form ( |2.13| ) treated in case 
B) above. Before proceeding further we remark that the case a = —2 in 
(p.40|) mentioned in Kamke [|][6.184, p 585] leads to p{w) = kw'^ exp{Ciw) 
= w and the resulting relation cannot be inverted for w = w{t). Equation 
(|2.40|) is solved and, as in ( p.l4|) , we obtain 



(2.40) 



Pi = I (1 - ^l-A{a + 2)d) , p2 = i (1 + ^l-4:{a + 2)d) , 4(a + 2)d ^ 1 

(2.41) 

Consequently due to (p.40|) 



w 



-pi/(a+2) 



pl 



-l/(a+2) 



dw = t + C. 



(2.42) 



If we set 



(p2-pi)(a + 2) 



a + 2 - pi 
the integral, /, in ( p.42| ) takes the form 

a + 2 /■ , ^7 l-l/(a+2) 



l-4(a + 2)c/ 



7(« + 3/2) 

a + 2 -7/2' 

(2.43) 



- / [Ci + C2']^'/^"+'Mt; = t + C 



a + 2- pi 
and, by virtue of the substitution 

Ci + C2V^ = 



(2.44) 



(2.45) 
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yields 



a + 2 



CsJl - 4(a + 2)c/ 



, Co 



r-l 



du = t + C, 1/7 



(2.46) 



The only possibility existing for / in ( ^^.46] ) to be computable is when both 
a and r are rational and in particular ||^ 

i) r = integer 

ii) (a + l)/2 = integer 

iii) (a + l)/2 + r = integer. 
1, we arrive at the known cases a = 



If r = 
case is a 



1 and r 
/ = 



2, ie 7 
3 



-2 or d = 0. The next simplest 
1/2. Then from ( p35|) , (|2:43| ) and (g^ 



121' 



12rf 
d 



u 



Co 



-2 + 25c/62 



1. 



C 



From (U), (|;34D, (|2^), (|2^ and (CT) with a 
deduce that 

20 



1 and 7 



u 



2u'Ci - —Ci{t + C) = Q 



33 



= - 



3m2m 



(2.47) 

1/2 we 
(2.48) 
(2.49) 



Relations (|2.48|) and ( p.49|) constitute an exact solution of ( p.33| ) provided 
the conditions ( p.47|) are satisfied. Equation ( p. 48 ) can be solved for u{t) and 
along the same lines as in Case B) above we arrive at a singularity pattern 
similar to that of (|2.22| ), ie branch points and essential singularities. 

We think that there is no essence in pursuing this further by examining 
other cases in i) - iii) above since we have already obtained a very general 
singularity behaviour of the solutions to ( |2.33 ) through ( 2.48 ) and ( 2.49| ). 

Now, in contrast to case B) above, some further exact results can be 
extracted for ( p.33| ) by following Kamke p [loc cit] and we obtain for ( p.33| ) 



P 



X, 



P = p{x) = XUi, Ui=Ui{ti), ti 

u'lUi + ul{a + 2) + bui + c = 



In a; 



(2.50) 
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whereas, after a straightforward calculation, 

b c 



a + 2' 



b' 



a + 2 



log {ul + a'ui + 6')^^^ — log 
= log \{Ci/x 



a+2' 



Ci = cons tant. 



27[2Va'^-46'] 



(2.51) 

For ( p. 511) to be of any use at all the quantity a' /[2^/a/^'^^4]y] has to be 
an integer and we consider in the following the simple case 



1^/ 

2 



a = V a'2 - 46' 



(2.52) 



and, therefore, (|2.51|) finally yields 



— uw + \a' = 

I 3 / I I \2a+4 

M = Ml + -.a , W = [Ci/x) 



(2.53) 



Along the lines of case B) above and by applying (|2.5CI| ) we deduce 



/ 2 
X = X < 



3a _j_ 
4 



4 ^ 



CI 



6a + 12 



+ 



1/2- 



9^6q+12 I ig 

^6a + 12 , '\2 \ ^/'^ 

9j.6a + 12 -r -^g 



1/3 ^ 



1/3 



(2.54) 



On integrating ( ^.54] ) we obtain 



3a _|_ 

4 



-4 + + ^ 



16 



1/2- 



1/3 ^ 



+ 



a 
' 4 



/^6a+12 ^ (oi 



16 



1/2- 



1/3 



> dz = t + C2 



(2.55) 



The integral in (|2.55| ) can be brought into a form more amenable to 
treatment by setting 



6a+12 



+ 



^ X 1/2- 

16 . 



1 1/3 






a 1 


+ 


"4 " i 



6a+12 



+ 



16 J 



1/3 



= y 

(2.56) 
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whence 



2a + 3 4a + 7 
' 2a+4 7/ 2a+4 ( 



-3) 



1 

'2a + 4 



{-y^-^)2a+4y 



2a + 5 - 
' 2a+4 



(2.57) 

If we reflect a little on (|2.57|) we conclude that only for specific a-values 
can we attempt to carry out the integration in ( p.57| ). For example, if a = 
-3/2, then ( p3TD , (|3|) and (|235|) - ( ^371) yield 



f + 4 + log 
x(t) = 



fa-3b/2)''+->/3^ 



y 



4/96 



= kit + C2) 



(2.58) 



as an exact solution to (2.33), where C2 and k are constants, provided 



-3/2, 



36V8, & > 0. 



(2.59) 



We remark that if we replace ( 2.52| ) by |a' = —^/a'"^ — 46' we again arrive at 
(|2.58|) - ( |2.59D , the only difference being that ( p.59| ) also holds for 6 < 0. From 
( |2.58| ) it also follows that x{t) has a singularity at t = tg, being calculated 
from ( p.58| ) by setting y = (—6)^/^. 

Evidently we can find other a-values, for instance a = —11/6, for which 
the integration in ( 2.57|) is possible and we would arrive at formula similar 
to ( [2.581 ). Likewise we may take into account other possibilities in place of 
(|2.52|) . Nevertheless there is little point in carrying this investigation further 
since it is almost certain that the resulting relations equivalent to ( [^.58D will 
not be invertible for the auxilliary function y{t). 



dy 

y-^ 



3 Painleve analysis of the second order 
equation 



We determine the leading order behaviour of ( 1.18| ) by means of the substi- 
tution X = ar^ and find that the value of p is —1 in accordance with the 
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coefficients of the rescaling symmetry, G2, and that the coefficient, a, is a 
solution of the equation 

a^- 3a + 2 = 0, (3.1) 

ie a = 1,2. To find the resonances we make the substitution x = ax~^ + 
f3x^~^. The condition that j3 be arbitrary is 

+ 3{a- l)r + 3a-4 



0, 



(3.2) 



ie r = —1,4 — 3a. With the two possible values of a we have r = —1,1 
for a = 1 and r = — 1, —2 for a = 2. The former gives the standard Right 
Painleve Series, representing a Laurent series in the neighbourhood of the 
moveable singularity, and the latter a Left Painleve Series, representing a 
Laurent expansion valid over the exterior of a disc surrounding the moveable 
singularity 

In terms oi t = t — the explicit expressions for the series are 



1 



x{t) = - [l + aiT - afr^ + afr^ - afr^ + . . . 



(3.3) 



for the Right Painleve Series and 

x{t) = 2r^ + ait''^ + a2t~^ - \ai (al - 302) t"^ + ^2)^^"^ (3.4) 

for the Left Painleve Series. For the latter we have written the expansion in 
terms of t instead of r because of the asymptotic nature of the expansion [^]. 

We can reconcile the two Laurent series, ( p.3|) and (|3.4| ), with the exphcit 
solution, (O) with 6=1, given in the previous section. In the case of the 
Right Painleve Series we have 



x{t) 



1 


2 - 


r 




1 


2 - 


r 





1 — OiT 
1 



1 + aiT. 
1 + 2aiT 

r(l + 2air) 



2 2 



.3^3 



4 4 I 



(3.5) 



which is just (|1.9| ) written slightly differently. If we now expand this as 



x{t) 



to 
t 



ai 
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2t- 



+ 



2tn- 



ai 



(--2 



2ti 



3 

^''0 ' 2 3 

, Oil '^1 'J^l y 



ai 



^-3 



(3.6) 



we can verify that the expansion is identical to that in ( [4 .111) when we rename 
the coefficients of and in (|3.6|). 



4 Painleve analysis of the two-dimensional 



system 



To determine the leading order behaviour of the system, ( |1.19D and (|1.2CI|) , 
we make the substitutions 



X 



(4.1) 



to obtain 



apr'' 



-(1 + a)(l + 2a)aV2p - 3(a + l)a/3r^'+^ - /^V^^ 



so that the powers to be compared are 



-1 p q 
1 2p p + q 2q 



(4.2) 



(4.3) 



in which we have removed a common p from the first line. If the first term in 
the right side of ( [4. 21 ) were zero, we would also be able to subtract a common 
q from the second line. Evidently there are three possible systems which 
should be analysed separately. 



4.1 The general case 

In general the product (1 + a)(l + 2a) ^ 0. In this case we can have two 
possible types of leading order behaviour. The first is the generic one for 
which p = —1 = q. The second is g = —l,p > —1. 
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We commence our analysis with the first case. The two coefficients a and 
/? are found from the solution of 



— a = ac? + a/3 
-/5 = -(1 + a)(l + 2a)a2 



3(1 + a)a/5 - /J^ 



(4.4) 



and we obtain a = 1,2 and /3 = — (1 + a),— (1 + 2a) respectively. 
To determine the resonances we make the substitution 



X 

y 



r-l 



VT 



(4.5) 



which allows both sets of coefficients to be treated at once. The condition 
that the coefficients /i and v be arbitrary is given by 



r — aa 
[1 + a){aa + 2a 



3) r-(3 



-a 

3a 



aa 







(4.6) 



with the solution r = —1,4 — 3a so that in the case a = 1 the resonances 
are at r = —1, 1, indicating a Right Painleve series, and in the case a = 2 
r = — 1, —2, indicating a Left Painleve series. Note that the resonances are 
independent of the value of a. The coefficients at the resonances are given 
by 



0", 



(4.7) 



where a is a parameter, for the cases a = 1 and a = 2 respectively. 

The other possible leading order behaviour is given by g = — 1 and p > — 1 
which means that the series for x starts at the constant term and so the 
standard Painleve analysis cannot be used. We write 



X 



aiT 



y 



biT 



(4.8) 



in which the summation convention is adopted and the counter i commences 
at zero. We substitute ( |4.8| ) into the system ( |1.19| ) and ( |1.20D and obtain 



lajT 



i-l 



aaiajT 



l)hT 



i~2 



-(1 + a)(l + 2a)aiajT 



(4.9) 

3(1 + a)a^bjT'+'-^ - hh,T'+^-^. 

(4.10) 
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We equate the coefficients of like powers of r to zero and obtain 



ao = 
ai = Oi 
02 = 

«3 = 1(3 + 4a)af 
04 = 



bo = l 

61 = 

62 = 3(1 + a)ai 

h = o 

64 = -i(l + a)(19 + 20a)a2. 



(4.11) 



It is a simple matter to verify that the ffist few terms of the series for x give 
the ffist few terms of the series for y when substituted into ( |1.21|) . 



4.2 The case (1 + a)(l + 2a) = 

As there are two roots to this quadratic equation, we have to consider two 
subcases. 



4.2.1 l + 2a = 

The two-dimensional system, ( p..l9| ) and ( |1.20| ), is now 

X = —^x^ + xy 

y = -Ixy-y^. (4.12) 

For the case when all terms are dominant, viz p = q = —1, the coefficients of 
the leading order terms are found from the solution of the system 

— a = —\o? + a(i 

-/5 = - (4.13) 

and are a = 1 and /3 = — ^. The resonances are given by r = ±1 and so we 
obtain a Right Painleve series. The arbitrary vector at the resonance r = 1 



o. (4.14) 



In addition to all terms contributing to the dominant behaviour there is 
the possibility that just some of the terms can be dominant. Such terms are 
termed subdominant. These terms must have the same symmetry which will 
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be more restrictive than the rescahng symmetry of the system as a whole 
and this property is termed subselfsimilarity. The first possible subdominant 
behaviour is given byp=— l,g>— 1. As the series for y does not contain a 
singularity, we must, as we did above, substitute a full series into the system 
( [4.12|) . We assume that 



X = aiT' ^ y = biT' (4.15) 

with the same conventions as before. We obtain the coefficients 

ao = 2 bo = 

ai = bi = 

02 = 62 = (4.16) 

03 = 63 = 

04 = 64 = 

so that we only obtain an isolated solution. 

There is another possibility and that is that the subdominant behaviour 
indicates the existence of an asymptotic series, ie we have p = —l,q < — 1. 
We write 

X = a,t-'-^ y = bit-'-^ (4.17) 

with the counter beginning at zero. Note that we do not expand in terms of r 
since the expansion is asymptotic. We obtain the following set of coefficients 

ao = 2 &o = 

''^1^ 0. ^^3^ . (4-18) 

0-2 = - 26i O2 = f oiOi 

a^ = \{a\~Sbi)ai 63 = | (a? - 61) 61 



and, if the first few terms of the series for x are substituted into ( |1.21|) with 
th 

y- 



this value of a, we obtain verification of the first few terms of the series for 



4.2.2 l + a = 

The system is now 



X = —x^ + xy 
.2 



-y- (4.19) 
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In the case that all terms are dominant the coefficients of the leading order 
terms are given by 

— a — —o? + cv/S 

-/3 = -0^ (4.20) 

and we find that a — 2 and (3 = 1. The resonances occur at r = — 1, — 2 and 
consequently we have a Left Painleve series. 

One of the possible patterns of subdominant behaviour is when q = 
— l,p > —1. We write 

x = aiT' y^hr'-^ (4.21) 
with the usual convention and find that we obtain the coefficients 



ao 


= 


bo 


= 1 


Ol 


= 


bi 


= 


02 


= 


b2 


= 


03 


= 


b3 


= b. 


04 


= 


&4 


= 


05 


= 




= 


Qq 


= 


be 


- 



(4.22) 



which is very strange since we are generating a solution for y although it is 
defined in terms of a zero function. 

If we try the pattern q = —l,p < —1, we find that the first few coefficients 

are 



(4.23) 



from which it is evident that we have the solution to the system to be 

X = 

1 hi hi hi 

f = t + p + ? + ? + - 

^ (4.24) 



ao 


= 


&0 


= 1 


Ol 


= 


bi 


= 61 


02 


= 


b2 


-bl 


03 


= 


b?> 


= hl 


04 


= 


&4 


= b\ 



t-hi 

which is a one parameter solution of more than slightly suspect appearance 
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5 Comment 



We can make several observations of the results of these calculations. 

1. There is a close relationship between the possession of the Painleve 
Property and the ease of explicit integration of this class of equations. 
Since they are all second order equations with two symmetries, they are 
reducible to quadratures and so integrable in the sense of Lie. However, 
as we have seen in Section 2, integrability in the sense of Painleve is 
a much rarer occurrence. We recall that integrability in the sense of 
Painleve imposes stronger requirements upon the solution since the 
function must be meromorphic. 

2. The properties of the two-dimensional system are independent of the 
value of the parameter a in the definition of the new variable except 
for the two particular values of a = — 1, — |. It is interesting to note 
that the quadratic equation which determines these values of a is the 
adjoint of the equation which determines the value of the coefficient, 
a, of the leading order term in the analysis of the original second order 
equation, (PTT|). 

3. In addition to the possibility of the existence of the Left Painleve Se- 
ries and the Right Painleve Series we have seen the possibility of the 
existence of an asymptotic expansion which really has nothing to do 
with the Painleve analysis even though it indicates singular behaviour. 

4. What is most disturbing is the existence of solutions which are appar- 
ently spurious. It remains an open question as to the correctness of the 
proper interpretation of these solutions. 
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